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We reconsider the Mott transition problem in the presence of long range Coulomb interactions. 
Using an extended DMFT, that sums an important class of diagrams absent in ordinary DMFT, we 
show that in the presence of Coulomb the Mott transition in two and three dimensions is discontin- 
uous as envisioned by Mott. 

PACS numbers: 71.30+h 
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The study of Metal Insulator transitions jjj (MIT) seen 
in a host of interacting electronic systems, has been one 
of the most challenging problems in solid state physics. 
These transitions can result from simple band structure 
or from electronic correlations. The former occurs even 
in nearly non-interacting systems if the conduction and 
valence bands split as some parameter like composition 
is varied. The correlation driven transition i.e., Mott 
transition, is however, the most interesting. The na- 
ture of these transitions depend on the interplay of band 
structure, magnetism and electronic correlations. The 
simplest scenario for a MIT was first put forth by Mott 
M . Mott suggested that a crystalline array of hydrogenic 
atoms i.e., atoms with just one electron in their outer- 
most shell, shows a zero temperature transition from a 
metal to an insulator as the density is decreased or effec- 
tively, as the distance between atoms increases Later, 
other possibilities for the MIT were put forth, namely, 
Slater's band transition Q where the system becomes 
insulating because of the doubling of the unit cell due 
to antiferromagnetism and the Brinkman-Rice scenario 
within the Hubbard model [|| which is characterised by 
a strong mass enhancement. 

Mott's original argument depended on the presence of 
coulomb long range interactions and pointed out that 
since the number of " free" electrons could only vary dis- 
continuously across the transition, the transition is nec- 
essarily first order. We recapitulate Mott's ideas below: 
the main premise was that, due to the long range nature 
of the interactions, the electrons and the holes in the two 
Hubbard bands will always form bound states and not 
exist independently. The condition for the formation of 
atleast one bound state, however, depends on the screen- 
ing of the coulomb interaction. Using Thomas-Fermi es- 
timates for screening by N electrons per unit volume, 
Mott showed that no bound states exist provided 



N 1/3 a H > 0.4 



(1) 



where an is the Bohr radius. This implies that the sys- 
tem is metallic as long as the density satisfies ([!]). For 
densities where ([J) is violated, the electrons and holes 
in the Hubbard bands always form bound pairs resulting 
in an insulating behavior. This implies that the number 
of carriers jumps at the transition. Since this results in 



a kink in the free energy, the transition is discontinuous 
and first order. For an illuminating discussion of these 
ideas see Ref. ||. 

In contrast, various studies of the Hubbard model 
HUJt]] indicate that the transition is continuous. The 
Hubbard model for interacting electrons has long been 
used as a prototype to describe various aspects of real 
systems jy . This model contains a term which describes 
the hopping of electrons between different atomic sites 
and another describing the coulomb repulsion felt by two 
electrons on the same atom. In the past few years, the 
nature of the MIT in this model has been clearly eluci- 
dated using the Dynamical Mean Field theory (DMFT) 
in the limit of infinite lattice coordination d Q]. It was 
shown that the MIT was characterised by a continuous 
vanishing of a Kondo like resonance at the Fermi level in 
the metal at the transition, leading to an insulator with 
a preformed gap. This led to the effective mass of the 
quasi-particles and hence the linear coefficient diverging 
from the metallic side at the transition HfjJ . Despite the 
discontinuous opening of the Mott Hubbard gap the de- 
struction of the metal at zero temperature treated with 
DMFT, is a continuous one. 

The purpose of this work is to incorporate some of the 
effects of the long range coulomb interactions and Mott's 
ideas into the framework of Dynamical Mean Field The- 
ory. For this purpose, we explore the effects of the long 
range interaction on the MIT seen in the single band 
Hubbard model j|, using a simple extension of DMFT 
iPUP^ . It was shown, in a model of spinless fermions, 
that this approach captures important 4 corrections [ pT| . 
This method was independently developed in Ref. |12| 1 
and it was applied to the problem of the breakdown of 
Fermi liquid theory. We first describe the approach by 
isolating a class of diagrams which can be formally con- 
trolled by scaling the interactions and the kinetic energy 
appropriately and which can be summed using impurity 
models. We then demonstrate that when this extended 
DMFT is applied directly to a 3 or 2 dimensional lattice 
with interactions having the coulomb form, it changes 
the Mott transition which was continuous in ordinary 
DMFT, to a discontinuous first order transition as envi- 
sioned by Mott. 
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The effective hamiltonian we use to describe our sys- 
tem is a generalization of the Hubbard hamiltonian 

H = tijCiaCj* + h - C - + V C t* C ^ 

(ij)a i<y 

+ y j Unj-\riii + ) Vij : nio-njv' : (2) 

i i^j^aa' 

The first term is the hopping matrix element for an 
electron from site i (representing the ion at Ri) to its 
neighboring site j, \i is the chemical potential, U is the 
coulomb repulsion felt by the electrons when they are on 
the same atom (i = j) and the normal ordered last term 
is the coulomb interaction between electrons on differ- 
ent atoms. The coupling constants are given by overlap 
integrals involving a set of chosen basis vectors like the 
Wannier or Hartree basis. For example, the on-site in- 
teraction U is given by 

U = J drdr'\u{r)\ 2 \u{r')\ 2 ^-^ (3) 

r e 2 
Vij = J drdr'\u{r - R t )\ 2 \u{r' - R^y—^ 

The hamiltonian (^|) is now studied in the dynamical 
mean field approximation. The first step in the DMFT 
is to scale the parameters in the large d limit appropri- 
ately such that the corresponding energy terms remain 
finite. In addition, the scaling should be chosen such that 
the terms of interest remain relevant in the large d limit. 
We adopt the following scaling: is scaled as \fd 3 , 

U — > U and Vij is scaled by Vd 1 3 . The leading dia- 
grams can then be summed up using the cavity method 
described in Ref. This scaling is a prescription for 
choosing a set of diagrams that contribute to the self en- 
ergy and hence the Green's function. DMFT prescribes 
a scheme for choosing different sets of diagrams that con- 
tribute to the different Green's functions. In the case of 
the Hubbard model (V = 0), this corresponds to retain- 
ing only skeleton diagrams constructed from U and the 
local Green's functions Gu in the self energy, resulting in 
the self energy £y being local. 




FIG. 1. Reduction of the self-energy diagrams in the limit 
of infinite d. 

In the presence of the longer range interaction, the above 
scaling retains all skeleton diagrams constructed using 
the local Greens function Gu and the interaction vertices 
U and Vij , such that every point i which has a vertex Vij 
originating from it has another vertex Vki terminating 



at it. In effect, this corresponds to replacing the U by 
an U e f f in the local self energy evaluated for a Hubbard 
model. An example of the diagrams retained is shown in 
Fig.|l|. We mention that the Hartree term (arising only 
from U) is generated in the cavity method. Depending 
on the basis chosen to derive the effective model param- 
eters (^) care should be exerted to see that the Hartree 
term is not double counted. We neglect the Fock term in 
the ensuing calculations , since the Fock term is of higher 
order in k than the Hartree term. 

Using the cavity method and integrating out all sites 
save a chosen site or cavity o, we obtain the following 
local effective impurity action with retarded interactions 

S eff = j drdr'^tA^GoHr-r'Mr') 

+J7n 0T (T)n 0i (r) - ^^(t)^ 1 ^ - t> 0<j <(t') (4) 

<J<j' 

where 

Oo 1 =d T -J2 t o t t oAf(r-r') (5) 

ij 

and the retarded interaction 

H^(T-T') = J2VoiVojIif(r-r') 

ij 

= J2V(q) 2 Tl^( q ,T-T') (6) 

9 

with G\j denoting the single particle Green's function 

and n^ 0) = (J2aa> n ^( T ) n 3<?' ( T '))o- The superscript (0) 
implies that the quantities have been evaluated in the sys- 
tem from which the site and all the links to this site have 
been removed. These in turn are related to quantities 
evaluated on the entire lattice. Summing over the sites 
leads to Go and no being determined by self-consistent 
equations involving local quantities evaluated on the full 
lattice. Note that the effect of all non-local quartic in- 
teractions is to dynamically screen the on-site repulsion, 
with the screening potential given by (§). These equa- 
tions depend on the nature of the lattice used. For ex- 
ample, on the Bethe lattice, g\j = Gij and = Hy. 
Retaining only the on site and nearest neighbor interac- 
tions, the self-consistency condition that should be sat- 
isfied by the retarded interaction on the Bethe lattice 
takes the simple form n^ 1 = V 2 Hi oc , where Hi oc is the 
local density density correlator. For arbitrary lattices 
and interactions (i.e. general tij and Vij ) the (extended) 
dynamical mean field equations can be derived by gener- 
alizing the discussion presented here, to the case where 

i — | ^ — 1 

the hoppings iy are scaled by v d 

To obtain the Green's functions of ([!]), we now define 
certain irreducible quantities. The lattice Green's func- 
tions, can be expressed in terms of a self energy £ which 
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is two particle irreducible and which becomes local in the 
limit of infinite dimensions. 



G(iuj n ,q) 



1 



E(iw„) 



(7) 



where e q is the dispersion on the lattice and iu) n are the 
Matsubara frequencies. Similarly, the density density 
correlator Tl(q, iu> n ) on the lattice defines an irreducible 
part II via the Dyson equation 



U(q,iUn) 



II 



l + V(q)U 



(8) 



where II is the sum over all polarization diagrams con- 
structed with the full G and interaction vertices V(q) 
such that all the diagrams are irreducible with respect 
to V(q). Within the infinite d approximation, since all 
vertex functions become independent of momenta, IT also 
becomes independent of momenta. The impurity model 
(EI) allows us to compute all local quantities, in particular 
£ and II and hence, the local density density correlator 
Tlioc and the local one particle Green's function Gi oc as 
functionals of the Weiss fields Q and H . 

n^{n O) 0o}(^n) - n-Hno.&KiwO - n^ l (iu n ) (9) 

and 

£{n o ,0o}(£w n ) = Q^i^n) - Gi oc {n ,g }(iuj n ) (io) 

Using (||) and (0), we can eliminate ft to obtain a self- 
consistent equation for II; oc = ^2 H(q,iu} n ) 

n ioe {n , g } = ^ 1 (11) 

Similarly, we obtain an equation for Gi oc 



Gi oc {Ho 7 Go} = 

* — * it, )_. — — 



„ -e q - £{n ,<5o} 



(12) 



These equations where derived by an appropriate scal- 
ing of the interactions and the hopping elements in the 
original model so as to obtain a well defined limit of large 
coordination. In the spirit of DMFT M, we can, however, 
regard these equations as defined on a finite dimensional 
lattice by replacing V by the usual coulomb interaction 
on the lattice. We use these equations to make qualita- 
tive predictions of the effect of the coulomb long range 
interaction in finite dimensions, on the order of the Mott 
transition seen. 

In earlier DMFT studies of the Hubbard model the 
continuous Mott transition was found to occur at a crit- 
ical U = U C 2 |r|] and was signaled by the vanishing of 
a Kondo like resonance at w = at U C 2- The continu- 
ous character of that transition was established using the 
projective self consistent method M. 



We briefly review the projective method below. Then, 
using the results of Ref. M, we study how long 
range Coulomb interactions treated within the extended 
DMFT, modify the previous discussion. The projective 
method uses the separation of two energy scales that ex- 
ists in the metallic phase close to the MIT i.e., wD where 
w is the weight of the Kondo resonance seen in the metal 
and scale of the Hubbard bands U. D is the half band- 
width. The high energy scales are then eliminated to ob- 
tain an effective theory governed by one low energy scale 
which is w, which goes to zero at the MIT Using this 
effective low energy model, which is a Kondo model of an 
impurity spin interacting with a bath of electrons, we can 
obtain the free energy or the ground state energy of the 
lattice problem close to the transition point. This free 
energy correctly describes the low energy and coherent 
part of the spectra. To order w, the low energy Kondo 
problem has been derived in Ref. M and is given by 



H = wTS.s L +Hb 



(13) 



Here Tib describes a band of low energy conduction elec- 
trons and Si represents the local spin operator of these 
electrons. S is the impurity spin and the Kondo cou- 
pling r is determined by matrix elements in the high 
energy sector (comprising the two Hubbard bands). The 
self-consistency conditions translate to conditions on the 
expectation value (S.s^). The high energy sector is an 
insulator and has a spin doublet ground state which is 
separated from the excited states by a large gap. Using 
the results of Ref. || , the simplest approximation to the 
high energy sector of the Hubbard model yields T — . 
Since we are interested only in the qualitative features of 
the transition, we replace the bath of electrons by a sin- 
gle electron. Therefore, using this toy model, the energy 
of the Kondo model is 



E K = --wT 



(14) 



Using the results of |TJJ , and taking into account the ki- 
netic energy of the lattice which is positive and of order 
w, the ground state energy of the lattice model has the 
following expansion in terms of w 

2D 2 

Eq = (a — (3 —jjr- )w + jw 2 + corrections (15) 

Here a, (3, 7 are all positive. The transition is determined 
by the U at which the coefficient of the term linear in 
uu vanishes. Since the energy minimum still occurs at 
w = 0, this transition is continuous. 

Though the DMFT and the self-consistent projective 
method are both formulated in infinite dimensions, it 
is nonetheless known that they the capture some as- 
pects of the physics of systems in finite dimensions rather 
well. We therefore, use the projective method to make 
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some predictions about the transition in three dimen- 
sions. First, in the presence of the long range coulomb 
interaction, there are corrections to (15), because of the 
screening of the on-site coulomb interaction given by 
(0) . This implies that U in ( p"s|) has to replaced by 
an U e ff. For the free energy (|15|), it is sufficient to 
consider the effective static on-site repulsion. Taking 
n(°) = IT + sub — leadingcorrections in (^) and using the 
fact that at zero frequency, II in (||) is by definition the 
compressibility k of the system, one obtains from (Q) 



U, 



eff 



u 



KV(q) 



(16) 



Note that V(q) is the Fourier transform of the coulomb 
potential in d dimensions. Since the effective interaction 
remains finite at the transition, the leading critical be- 
havior of the compressibility k is the same as that found 
in the Hubbard model and k is given by its value in the 
Hubbard model 



D 



(17) 



Note that K is finite in the metal and goes to zero 
smoothly at the transition. This implies that the effec- 
tive repulsion seen by the electrons in the metallic phase 
is reduced. In the insulator k = and the bare on-site 
repulsion is not screened at all. This is physically correct 
because one does expect the itinerant electrons in the 
metallic phase to better screen the coulomb repulsion. 
Using ( pj| ) in (16), and performing the q sum in d = 3 
we find that the effective repulsion felt on-site is 



U, 



eff 



U 



•,1/2 



(18) 



where a is some positive constant. The screening is pro- 
portional to log it; in two dimensions. Substituting JlS| ) 
in Jl5|), we see that the expansion of the energy in terms 
of the low energy weight w has non-analytic terms i.e., 



2D 2 . 2D 2 s 

E = (a- P——)w - afj—pw- 
U U z 



■ 7«; 



correction 



(19) 



Notice that the parameter w in the projective self con- 
sistent method, neatly embodies the notion of "number 
of free carriers" in Mott's original number. The non- 
analyticity in (^) arises purely from the self-consistency 
conditions ( fjl| ) when coulomb interactions are taken into 
account. We now see that at the transition, though the 
linear term still vanishes, the minimum of the energy 
which formerly was at w = is now shifted to a non- 
zero value because of this non-analytic term. This is 



a feature of a first order transition! Hence, we see that 
taking into account the long-range nature of the coulomb 
interaction does make the Mott transition first order in 
d = 3. Regarding the effects of coulomb interactions in 
two dimensions, we find that the screening results in a 
term proportional to w log w in the free energy ([l5]), again 
making the Mott transition discontinuous. Even-though 
the DMFT is too crude an approximation in d = 2 where 
spatial fluctuations are large, we argue that it still cap- 
tures some of the essential physics of strongly correlated 
systems. 

To conclude, we have studied the effects of long range 
coulomb interactions on the Mott transition using an ex- 
tended DMFT. We find that the coulomb interactions 
dynamically screen the effective on site interaction. We 
find that to lowest order, the screening is determined by 
the low energy scale w which is related to the width of the 
Kondo resonance at the Fermi level. Consequently, the 
screening is zero in the insulator. Extending our analysis 
to three dimensions, we find that the screening term is 
non-analytic in w resulting in non-analytic terms in the 
free energy. As a result the Mott transition which was 
continuous in the absence of these dynamical screening 
terms, now becomes a first order transition vindicating 
Mott's ideas. 
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